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Abstract. There is compelling evidence that black holes (BHs) in cluster centers vigorously interact 
with their surroundings, indicating that any realistic model of cluster formation needs to account for 
these processes. Here we use high-resolution cosmological simulations of a large cluster and group 
sample to study how BHs affect their host systems. We focus on two specific properties, the halo 
gas fraction and the X-ray luminosity-temperature scaling relation, both of which are notoriously 
difficult to reproduce in self-consistent hydrodynamical simulations. We show that BH feedback can 
solve both of these issues, bringing them in excellent agreement with observations, without alluding 
to the 'cooling only' solution that produces unphysically bright central galaxies. 
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INTRODUCTION 

Hydrodynamical cluster simulations which only include radiative cooling typically suf- 
fer from excessive overcooling within the densest cluster regions. This results in a very 
large fraction of cold gas and consequently a large amount of stars that would form 
out of it. While the associated removal of low-entropy gas breaks the self-similarity of 
the cluster scaling relations in a way that resembles observations [HiSIl], this 'cooling- 
only' scenario is generally discounted because observed clusters do not show such strong 
cooling flows and enormously bright central galaxies. Instead, it is widely believed that 
some non-gravitational energy input strongly affects the intracluster medium (ICM). 

There has been considerable effort [see 4, and references therein] to include feedback 
mechanisms associated with star formation in hydrodynamical simulations in order 
to reduce excessive overcooling in cluster cores. However, so far simulation models 
have failed to simultaneously reproduce the masses and colors of central galaxies, the 
observed temperature and metallicity profiles, as well as the observed X-ray luminosity- 
temperature (Lx — T) relation. In particular, the X-ray luminosities have been found to 
be substantially larger than the observed values ^ on the group scale. 

Feedback from active galactic nuclei (AGN) might resolve the discrepancies by pro- 
viding a heating mechanism that offsets cooling, lowers star formation rates and removes 
gas from the centers of poor clusters and groups, thereby reducing their X-ray luminosi- 
ties to values compatible with observations. We investigate whether AGN feedback can 
indeed solve these problems by performing simulations of galaxy clusters and groups 
that employ a state-of-the-art model [5] for BH growth and feedback. 



THE SIMULATIONS 



We have performed high-resolution cosmological hydrodynamical simulations of a large 
galaxy cluster and group sample, accounting for radiative cooling, UV heating, star 
formation and supernovae feedback. For each halo, two kinds of simulations were 
performed. One containing the physics just described and an additional one including 
a model for BH growth and associated feedback processes as in [6J and [5]. This allows 
us to compare the very same clusters simulated with and without AGN heating in order 
to clearly pin down the imprints of AGN activity on cluster and group properties. 

The simulations are described in detail in [7]. Here we just summarize the main 
features of our BH growth and feedback model. We assume that low-mass seed BHs are 
produced sufficiently frequently that any halo above a certain threshold mass contains 
one such BH at its center. In the simulations, an on-the-fly friends-of-friends group 
finder puts seed BHs with a mass of IO^H^^Mq into halos when they exceed a mass 
of 5 X lO^''/i^^M0 and do not contain any BH yet. The BHs are then allowed to grow by 
mergers with other BHs and by accretion of gas at the Bondi-Hoyle-Lyttleton rate, but 
with the Eddington limit additionally imposed. 

We use two distinct feedback models depending on the BH accretion rate (BHAR) 
[see 111]. For large accretion rates above 0.01 of the Eddington rate, AGN feedback is 
assumed to occur in a "quasar-mode", where only a small fraction of the luminosity is 
thermally coupled to the ICM. We adopt this thermal heating efficiency to be 0.5% of 
the rest mass-energy of the accreted gas, which reproduces the observed BH mass-bulge 
stellar velocity dispersion relation [8]. For lower BHARs, we assume that feedback is 
in a so-called "radio-mode", where AGN jets inflate hot, buoyantly rising bubbles in 
the surrounding ICM. The duty cycle of bubble injection, energy content of the bubbles 
as well as their initial size are determined from the BHAR. We assume the mechanical 
feedback efficiency provided by the bubbles to be 2% of the accreted rest mass-energy, 
which is in good agreement with observations of X-ray luminous elliptical galaxies 

It was shown in [St] that this model leads to a self -regulated BH growth and brings BH 
and stellar mass densities into broad agreement with observational constraints. 

Finally, we obtain realistic X-ray luminosities and spectroscopic temperatures for 
each simulated halo, as described in [0]. 

HALO GAS FRACTIONS 

In Fig. [H we show the ratio of gas mass to total mass inside the radius rsoo as a 
function of halo X-ray temperature. The most obvious effect of the AGN feedback is 
the significantly reduced gas fraction at the low temperature end of our sample, i.e. in 
poor clusters and groups. There the AGN heating drives a substantial fraction of the gas 
to radii outside of rsoo- This lowers halo gas fractions in spite of the reduced fraction 
of gas that is converted into stars in the runs with AGN. The potential wells of massive 
clusters are, on the other hand, too deep for AGN to efficiently remove gas from them. 
Thus the effect of the suppressed star formation becomes more important towards more 
massive systems. While the gas fraction in the very inner regions of massive clusters is 
somewhat reduced by the AGN, we find it unchanged or slightly increased within rsoo- 
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FIGURE 1. Halo gas mass fractions within rjoo of clusters and groups simulated without (circles) and 
with AGN feedback (squares). The arrows illustrate the effect of the AGN heating for each halo. For 
comparison, we show constraints obtained from X-ray observations. 



THE Lx-T RELATION 

By removing gas from the centers of poor clusters and groups, AGN heating also 
suppresses their X-ray luminosities and affects the Lx — T relation. In Fig. |2l we plot the 
X-ray luminosities L500 against the spectroscopic temperatures T^qq, for all our halos. 

Without AGN, we obtain substantially larger X-ray luminosities for poor clusters 
than observed, while for massive clusters there is reasonable agreement. This finding is 
in line with previous studies [4], and is a manifestation of the long-standing problem 
to explain cluster scaling relations in hydrodynamical simulations. However, when we 
include the AGN feedback the discrepancies between simulated and observed Lx — T 
relation are resolved. In particular, the relation is steepened on the group scale, as AGN 
heating removes a larger fraction of gas from smaller halos and thereby reduces their 
X-ray luminosities. For massive clusters, the effect of AGN is less important. Overall, 
the Lx — T relation obtained from the simulations with AGN feedback is consistent with 
observations at all mass scales, from massive clusters to small groups. 



SUMMARY AND CONCLUSIONS 

We have performed very high-resolution simulations of a galaxy cluster and group 
sample and compared simulations with and without AGN feedback. We find that: 

• AGN feedback significantly lowers the gas mass fractions in poor clusters and 
groups, even though fewer baryons are turned into stars. This is because the AGN 
heating drives gas from halo centers to their outskirts. Both the gas and stellar 
fractions are in a much better agreement with observations when AGN are included. 

• AGN feedback significantly reduces the X-ray luminosities of poor clusters and 
groups, while the X-ray temperature stays roughly the same or is even slightly 
reduced. This results in a steepening of the Lx — T relation on the group scale. 
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FIGURE 2. X-ray luminosity as a function of spectroscopic temperature for clusters and groups 
simulated without (circles) and with AGN feedback (squares). For each halo, the arrow illustrates the 
effect of the AGN heating. Data from observational X-ray studies is shown for comparison. For a subset 
also the best-fit power-law Lx — T relation is plotted. Including AGN feedback into the simulations 
drastically improves the agreement with observations and resolves the discrepancy on the group scale. 



• The Lx — T relation obtained from simulations with AGN feedback is in excellent 
agreement with observations at all mass scales. 

We find it extremely encouraging that this simple model for BH growth and feedback 
is capable of bringing the analyzed properties of galaxy clusters and groups into much 
better agreement with observations. This resolves a long-standing problem in their 
hydrodynamical modeling and opens up new exciting possibilities for using simulations 
to investigate their properties and their co-evolution with central supermassive BHs. 
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